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1—5 ' Abstract 

In this paper we develop a tractable framework for SINR analysis in downlink heterogeneous cellular 

networks (HCNs) with flexible cell association policies. The HCN is modeled as a multi-tier cellular network 

where each tier's base stations (BSs) are randomly located and have a particular transmit power, path loss 

Q , exponent, spatial density, and bias towards admitting mobile users. For example, as compared to macrocells, 

picocells would usually have lower transmit power, higher path loss exponent (lower antennas), higher spatial 

^ ■ density (many picocells per macrocell), and a positive bias so that macrocell users are actively encouraged to 

^— V ' use the more lightly loaded picocells. In the present paper we implicitly assume all base stations have full 

cn 

[■^^ . the whole network or a certain tier, which is equivalently the downlink SINR cumulative distribution function. 

The results are accurate for all SINRs, and their expressions admit quite simple closed-forms in some plausible 
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special cases. We also derive the average ergodic rate of the typical user, and the minimum average user 
throughput - the smallest value among the average user throughputs supported by one cell in each tier. We 
observe that neither the number of BSs or tiers changes the outage probability or average ergodic rate in an 
interference-limited full-loaded HCN with unbiased cell association (no biasing), and observe how biasing alters 
the various metrics. 

I. Introduction 

Heterogeneous cellular networks (HCNs) comprise a conventional cellular network overlaid with 
a diverse set of lower-power base stations (BSs) such as picocells [[II, femtocells [[3, [[3l, [[H, and 
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perhaps relay BSs |I3. Heterogeneity is expected to be a key feature of 4G cellular networks, and 
an essential means for providing higher end-user throughput [IH, [|3 as well as expanded indoor and 
cell-edge coverage. The tiers of BSs are ordered by transmit power with tier 1 having the highest 
power. Due to differences in deployment, they also in general will have differing path loss exponents 
and spatial density (e.g. the number of BSs per square kilometer). Finally, in order to provide relief to 
the macrocell network - which is and will continue to be the main bottleneck - lower tier base stations 
are expected to be designed to have a bias towards admitting users |l6l, since their smaller coverage 
area usually results in a lighter load. For example, as shown in Fig. [U a picocell may claim a user even 
though the macrocell signal is stronger to the user. The goal of this paper is to propose and develop 
a model and analytical framework that successfully characterizes the signal-to-noise-plus-interference 
ratio (SINR) - and its derivative metrics like outage/coverage and data rate - in such a HCN with 
arbitrary per-tier association biases. 

A. Motivation and Related Work 

The SINR statistics over a network are, unsurprisingly, largely determined by the locations of the base 
stations (BSs). These locations are usually unknown during the design of standards or even a specific 
system, and even if they are known they vary significantly from one city to the next. Since the main 
aspects of the system must work across a wide variety of base station deployments, the BS locations 
are usually abstracted to some baseline model. The most popular such abstraction is the hexagonal grid 
model, which has been used extensively both in industry |[8l, flU, ifTOl and academia [fTT]|. [[T2|. over 
at least the past three decades and it remains ubiquitous up to the present day. The SINR expressions 
resulting from such a model are complex and depend on multiple random variables, and are so metrics 
of interest are usually estimated by Monte Carlo methods. This complicates the understanding of 
how the various system parameters affect such metrics. Furthermore, given the previously mentioned 
deployment trends, the continued relevance of the hexagonal grid model is open to debate. 

An alternative is to model the BS locations as random and drawn from a spatial stochastic process, 
such as the Poison point process (PPP). Such an approach has been advocated as early as in 1997 [fT3l . 
lfT4l . [fTSl even for traditional coverage-focused deployments. Recent work has found the complete 
downlink SINR distribution under fairly general assumptions, and shown that at least in one fairly 
regular macrocell deployment, that it is about as accurate as the grid model fT6l|. In the interim, the 
PPP model has been used extensively for modeling unplanned networks such as femtocells [[TtII . [[TSl . 
lfT9]l or ad hoc networks ll20l, EH; and for both it is now fairly well-accepted [l22ll. [l23ll. 



In the context of the heterogeneous deployments of interest in the present paper, a tractable model 
for SINR in a general K-ticr downlink HCN has recently been found under certain assumptions 
[[24|. Il25l . The results are exact for positive SINR (in dB) only, do not include biasing, and the per- 
tier coverage probabilities are not given. Clearly, it is desirable to develop a more general analytical 
model which is accurate for all SINRs allows for a flexible cell association. This would allow for the 
consideration of biasing in the overall system design, which is considered a very important factor in 
the overall performance of an HCN |l26l, [|27l, (H, [lH. We also note that independently, 1161 has 
been extended to the case of an unbiased two-tier network, where both macro BSs and femto BSs are 
located according to a PPP [301 . 

B. Contributions 

We present a novel analytical model and set of results for SINR in downlink HCNs with flexible 
cell association. The model characterizes the K tiers of a cellular network by the transmit power, 
BS spatial density, path loss exponent, and bias factor. We assume that in addition to path loss that 
there is Rayleigh fadingl|. Unlike ll24l . [|25l where each mobile user connects to the BS offering the 
highest instantaneous SINR, we assume the user connects to the BS that offers the maximum long- 
term averaged received power (with biasing), i.e. fading is ignored. This results in a completely new 
approach to modeling and analysis of the HCNs. Specifically, we provide the following theoretical 
contributions. 

First, we derive a pair of prerequisite quantities that are of interest in their own right. These are the 
(i) per-tier association probability, which is the probability that a typical user is associated with each 
tier, , and (ii) the distance to the serving BS for a given bias. These are both probabilistic functions. 
From the per-tier association probability, the average number of users associated with a BS in each 
tier is derived, which captures the effect of cell association on the cell load of each tier. 

Next, we derive the complete outage probability over all SINRs, with arbitrary biasing. The outage 
probability relative to an SINR threshold is equivalently the CDF of SINR for a randomly selected 
mobile in the network (or a certain tier). Although the results are not closed- form for the most general 
case, they are easily computable. In addition, they reduce to closed-form in several cases, for example 
an interference-limited network with the same per-tier path loss exponents and no biasing. In that 

'Lognormal shadowing on both the desired and interference signals was considered in 1161 for one- tier, but significantly degrades 
tractability while not changing the main trends. We neglect it here for expediency, but it could be handled in a straightforward manner 
following the approach of 1161 . 



particular case, the outage probability does not depend on the spatial density, transmit power or even 
the number of tiers. Intuitively, adding more infrastructure of any type has a negligible affect on the 
coverage or outage, so primarily serves to provide more throughput (since more users can be served 
simultaneously as the number of base stations increases). 

Finally, we derive two measures of spectral efficiency in our proposed HCN model: the average 
ergodic rate of a randomly chosen user in the whole network (or a certain tier), and the minimum 
average user throughput - the smallest value among the average user throughputs supported by one 
cell in each tier. The average ergodic rate, like the outage probability, is independent of BS transmit 
power, BS density, and the number of tiers K in interference-limited HCN without biasing, whereas 
it is strictly worsened by biasing in a fully-loaded (full queue) network. However, rather than drawing 
conclusions about the relative merits of different bias values, the main contribution of this paper is to 
provide a tractable SINR framework for doing so in future work, since it seems clear that the loading 
on each tier will significantly affect the cell association policies. 

II. Downlink System Model 

A fairly general model of an HCN would include K tiers of BSs that are distinguished by their 
spatial densities, transmit powers, path loss exponents, and biasing factors. For example, a high-power 
macrocell BS network is overlaid with successively denser and lower power picocells (with positive 
biasing) and femtocells (also with positive biasing) in a plausible K = ?> scenario as seen in Fig. [21 
Femto BS locations will generally be unplanned and so are well-modeled by a spatial random process 
[fTSl . [fT9l , [|30l . II3TII . Perhaps more surprisingly, macrocell BSs may also be reasonably well-modeled 
by a random spatial point process, with about the same as accuracy as the typical grid model [fT6l . flHl, 
[fT3l . Under this model, the positions of BSs in the jth tier are modeled according to a homogeneous 
PPP $j with intensity \j and users also are located according to a homogeneous PPP ^^^^ with intensity 
A*^"^ that is independent of {^j}j=i,-,k- 

The downlink desired and interference signals experience path loss, and for each tier we allow 
different path loss exponents {aj}j=i^...^K > 2, which is permissible in consideration with empirical 
value in cellular network ll32l . Rayleigh fading with unit average power models the random channel 
fluctuations. No intra-cell interference is considered, i.e. orthogonal multiple access is employed within 
a cell. The noise is additive and has constant power W but no specific distribution is assumed. Every 
BSs in the jth tier uses the same transmit power {Pj}j=i ... x- 

Denote k E {1, ■ ■ ■ , K} as the index of the tier with which a typical user associated. To evaluate 



the outage probability and achievable rate, we shift all point process so that a typical user (receiver) 
lies at the origin. Despite the shift, the distribution of potential interfering BS is still a homogeneous 
PPP with the same intensity. We thus denote |Yfcj| as the distance from BS i G $fc to the origin, i.e. 
the typical user. Also denote {-Rj}j=i ... x as a distance of the typical user from the nearest BS in the 
jth tier, which is a main parameter for the proof following Lemma 1 and 3. 

A. Flexible Cell Association and Cell Load 

We assume open access, which means a user is allowed to access any tier's BSs. This may not be 
true in some cases (notably femtocells) but provides best-case coverage, and we leave extension to 
closed-access to future work. We consider a cell association based on maximum biased-received-power 
(BRP) (termed biased association), where a mobile user is associated with the strongest BS in terms 
of BRP at the user. The BRP {P.^j^^i ... ^ is 

Pr,, = PMRj/roy^^B,, (1) 

where Lq is the path loss at a reference distance tq (typically about (47r/z/)~^ for ro = 1, where u 
denotes the wavelength). All BSs in the jth tier adopt identical bias factor Bj, which is positive value. 
When {Bj}j=i...K = 1^ i-C- no biasing, biased association is the same as conventional cell association 
(termed unbiased association), where the user connects to the BS that offers the strongest average 
power to the user. In particular, when Bj = l/Pj for all j = 1, ■ ■ ■ , K, the user connects to the BS 
with the lowest path loss. Employing Bj > 1 extends the cell range (or coverage) of the jth tier BSs 
as shown in Fig. [2l This allows the tier selection to be tuned for cell load balancing or other purposes. 
For clarity of exposition, we define 

p.A^§.AEia=^ (2) 

-Oc -Dfc Ofc 

which respectively characterizes transmit power ratio, bias power ratio and path loss exponent ratio of 
interfering to serving BS. Note that Pk = Bk = ak = I. 

Each tier's BS density and transmit power as well as cell association determine the probability that a 
typical user is associated with a tier. The following lemma provides the per-tier association probability, 
which is essential for deriving the main results in the sequel. 

Lemma 1. The probability that a typical user is associated with the kth tier is 



If {aj} = a, the probability is simplified to 

■^' -y^ x,{PAr"' ) ■ ^*^ 

Proof: See Appendix lAl ■ 

Lemma 1 confirms the common sense that a user prefers to connect to a tier with higher BS density, 
transmit power, and bias. We further observe that BS density is more dominant in determining Ak than 
BS transmit power or bias factor because - < 1 (a > 2) both in practice and in our network model. 
From Lemma 1, we easily derive the average number of users per BS in the kth tier, which quantifies 
the cell load of each tier. 

Lemma 2. The average number of users associated with a BS in the kth tier (also the cell load of the 
kth tier) is given as 

Mk = 27rA(^^ f r exp I -vr J2f=i XAh^if'"'^'^'^' } dr. (5) 

7/' {aj} = a, (O simplifies to 

A(") 

Mk = T7 ^^r^ . (6) 

A. + Ef=i,,y.A,(P,5,)2/" 

Proof: Denoting 5*, N, N^. ' , and NJ^' as the area of the entire network, number of users in the 
entire network, average number of kth tier BSs and kth tier users, respectively, we obtain NJ^' = 
AkN = AkX^^^S and A^^ = XkS. From the relations, the number of users per BS in the A;th tier is 
given by 

A4 = % = ^. (7) 

Combining © with ^ gives the desired result. In case of {aj} = a, the desired result follows by 
using (l39l) instead of ©. ■ 

Higher BS transmit power or biasing lead more users to connect with each BS of the corresponding 
tier. Comparing Lemma 1 with Lemma 2 admits the perhaps surprising observation that as the BS 
density of the A;th tier increases, more users are associated with the A;th tier but the number of users 
per BS decreases even in that tier. The decrease in the number of users per BS is because in (|7]), Ak 
scales like y^ less than A^ for a positive value of c (see (|39|)). On the other hand, as more users are 
associated with the fcth tier, less users are associated with other tiers. It implies that deploying more 
BSs reduces the cell load of the corresponding tier as well as other tiers, whereas employing higher 



transmit power or bias factor increases the cell load of the corresponding tier while decreasing that of 
other tiers. 

B. Statistical Distance to Serving Base Station 

We consider a typical user at the origin associated with the A;th tier. Denote Xk as the distance 
between the user and its serving BS. Since BSs are deployed as a PPP, Xk is a random variable 
described by its probability density function (PDF) given in Lemma 3. 

Lemma 3. The PDF fx^ (x) of the distance Xk between a typical user and its serving BS is 

Proof: Since the event of X^ > a; is the event of Rk > x, given the typical user's association 
with the kth tier, the probability of Xk > x can be given as 



nXk >x]= F[Rk >x\n = k] = ^"'^pj^!:^j ^^ , (9) 

where P[n = k] = Ak follows from Lemma 1, and the joint probability of Rk > x and n = A; is 



F[Rk > x,n = k]=F 



Rk > x,P,^k{Rk) > maxPrj 



(a) 



W, 



J X 

/■oo 

J X 

'2n\k f rexp(-nJ2f=i>^j{PjBjf/"^r^/"'\dr, (10) 

J X 

where (a) follows from ([U), and (b) is given from (|35l) and (|36l) . Plugging (fTOl) into ^ gives 

F[Xk >x] = — -^ / r exp { -n J2f=i XjiPjEj^/^^r^/^^ } dr, (11) 

The CDF of Xk is Fx^{x) = 1 - ¥[Xk > x] and the PDF is given as 

fxdx) = ^^ = ^xexp{-vrEf.,A,(P,P,)^/"^x^/""^} ■ 

■ 
Given a typical user associated with the fcth tier, the interference at the user at distance Xk = x is 
the sum of aggregate interference from all BSs, which are outside a disc of radius (P,Pj)^/"^x"j for all 
j = 1, ■ ■ ■ , K (see (a) in (flOl)). Note that this gives a large difference in the interference model between 
a cellular network with open access and previous studies for wireless ad hoc networks and two-tier 
femtocell networks with closed access, which had no such lower limit of the minimum distance from 
interfering BSs. 



III. Outage Probability 

Define the outage probability O as tlie probability that the instantaneous SINR of a randomly located 
user is less than a target SINR. Since the typical user is associated with at most one tier, from the law 
of total probability, the probability is given as 

K 

= Y,OkAk, (12) 

k=l 

where Ak is the per-tier association probability given in Lemma 1 and Ok is the outage probability of 
a typical user associated with kth tier. For a target SINR r and a typical user SINRfc(x) at a distance 
X from its associated BS, the outage probability is 

Ok = E,[F[Simk{x)<r]] (13) 

The metric is the outage probability averaged over cell coverage (defined by fx^i^) ^^ Lemma 2). It 
represents the average fraction of the cell area that is in outage at any time. The metric is also exactly 
the CDF of SINR over the entire network. 

The SINR of a typical user at a random distance x from its associated BS in tier k is 

SINRfe(x) = -^ Pk9k,oX-''' 

where gkfl is the exponentially distributed channel power with unit mean from the serving BS, \Yji\ is 
the distance from BS i G ^j\Bko (except the serving BS Bko) to the origin, and hji is the exponentially 
distributed channel power with unit mean from the ith interfering BS in the jth tier. 

A. General Case and Main Result 

We now provide our most general result for the network outage probability from which all other 
results in this section follow. 

Theorem 1. The outage probability of a typical user associated with the kth tier is 

o \ poo , . 

0. = l--^y^ xexp{-^-7rEf=iC^,a;2/s.}da;, (15) 

where SNR = ^''^^~"" and 



C, = \,p^'-^ 



Bp+Z{T,a,,B, 



with Z{t, aj, Bj) = J_^ 2-^1 



1 1_A;2-A;_^ 



, and 2-Fi [■] denotes the Gauss hypergeometric 



function. It follows that the outage probability of a randomly chosen user is 



= 1 -Y, 2vrA,. / x exp { -^ - TT Ef=i C,x^/^^ ] dx. 
k=i "^0 ^ ^ 



(16) 



Proof: See Appendix O ■ 

Although Theorem 1 does not give a closed-from expression, the integral is fairly easy to compute. 
For some plausible cases, we obtain simplified results in the following section. 



B. Special Case: Interference-Limited Network 

Since the BS density is typically quite high in heterogeneous cellular networks, the interference 
power easily dominates thermal noise. Thermal noise can often therefore be neglected, i.e. W = 0, as 
we do in the rest of this section. 

I) No Noise, Equal Path Loss Exponents {«j} = a: 

Corollary 1. When {oj} = a, the outage probability of the kth tier and overall network is respectively, 



OkiT,{\j},a,{B,}) = l 



T.U^APAf'" 



^K 



52/a 



2/q 



YJ^i^Pr B;"^ + Z{r,a,B,) 



(17) 



and 



K ( K 



A 



k=i lj=i 



Afe 



5f" + Z(r,aJ,) 



(18) 



Proof: If {«j} = a, then aj = 1. From (fT5l) we obtain 



Ok{T,{\j},a,{B,}) = l 



2ti\i 



„2 v^X 



2/« 



X exp <^ -vrx X]j=i Aj-P,- B-'" + Z(r, a, Bj) 



J2/a 



dx. 



Employing a change of variables x"^ = v and /q°° e~^^df = ^, and plugging (l39l) give the result in 
(fTTl) . Similarly ((Ml) leads to the result in ^^. ■ 

This expression is simple and practically closed-form, requiring only the computation (or lookup) 
of a Z[a, 6, c) value. When {Bj} = 1, i.e. no biasing, the probabilities are further simplified. 

Corollary 2. When {aj} = a and {Bj} = 1, i.e. unbiased association, the outage probability of kth 
tier and the network outage probability are given as 



1 



0(r, {A,}, a, 1) = Okir, {A,}, a, 1) = 1 - —— - 

l + Z[T,a,l 



(19) 
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Proof: Plugging Bj = 1 in (flTl) . we obtain 

Ofc(r,{A,},«,l) = l- 



1 + Z(r,a,l)' 

which is independent on the index k. From (fT2l) and Yl!k=i'^k = 1, we obtain O = Ok^^^-^Ak = Ot- 
This gives the result in (fT9l) . ■ 

This is much simpler than Corollary 1, in particular the outage probability is now independent of the 
BS transmit power Pj and BS density \j for unbiased association. It is also independent of the number 
of tiers K. This means that an open access unbiased HCN can randomly add new infrastructure, and the 
outage probability is not increased. This perhaps counter-intuitive result makes sense when one realizes 
that a mobile user always associates with the strongest BS, and so the SINR statistics do not change 
as more BSs are added. Interference management appears to be required to enhance downlink quality, 
since heterogeneous deployments do not change the outage probability. This "invariance" property 
is also observed in a interference-limited HCN with a different cell association ll25l . where the user 
connects to the BS offering the maximum instantaneous SINR, whereas here the user connects to the 
BS offering the maximum average received power. 

2) No Noise, Equal Path Loss Exponents {dj} = 4.- When {aj} = 4, the Gauss hypergeometric 
function in Corollary 1 and 2 collapses to a simple arctangent function denoted by arctan(). This 
provides a yet simpler expression for the outage probability. 

Corollary 3. If {aj} = 4, the outage probability of kth tier and overall network is respectively, 



Ok{T, {A,}, 4, {S,}) = 1 ^ ^^_ (20) 

Ef=i^j\JPjBA^ + Vr/fi,arctan(^r/i?,)) 
and 



0(r,{A,},4,{%) = l-f;|f^^^^(l + y^ 

fc=i Lj=i ^ ^ 




Proof: When a = 4, from (I44|) 
Z{t,A,Bj) = -y/r J°/g— YT^dw = i/r arctan(-\/ r/i?j). Combining with (flTl) and (fTSi) gives the 
desired results. ■ 

When {Bj} = 1, i.e. no biasing, the outage probability is simplified to a single trigonometric function 
and a single variable. 
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Corollary 4. When {aj} = 4 and {Bj} = 1, i.e. unbiased association, the outage probability of kth 
tier and the network outage probability are given as 

0(r, {A,}. 4. 1) = 0,(r, {A,}. 4. 1) = 1 - ^^^J,,,,,,^ . (22) 

Proof: When {aj} = 4 and {Bj} = 1, from dSl), Z{t, 4, 1) = ^/t J^^ r^d-u = a/t arctan(A/r). 
Combining with (fT9l) gives the desired results. ■ 

As the simplified expression of Corollary 2, this corollary also shows that the outage probability of 
each tier is the same for all tiers, and even it is the same as the outage probability of overall network. 
This implies that adding small pico and femto BSs to the macrocell network does not change the SINR 
distribution of each tier, because the increase in interference power is counter-balanced by the increase 
in signal power. We thus expect that both network sum throughput and per-tier sum throughput linearly 
increases with the number of BSs without any interference management techniques, since the SINR of 
each tier cell remains same. Furthermore, for a single-tier K=l, the expression of (|22l) is the same as 
the result of 11161 eq. (25)]. Our work thus extends the single-tier approach [I16B to a general multi-tier 
in the special case of Rayleigh fading. 

IV. Spectral Efficiency 

We derive the average ergodic rate and the minimum average user throughput to measure spectral 
efficiency performance of the network. The average ergodic rate is obtained with similar tools as 
the outage probability was in Section III. The minimum average user throughput is derived from the 
average ergodic rate. Both metrics are computed in units of nats/sec/Hz to simplify the expressions 
and analysis, with 1 nat/s = 1.443 bps. 

A. Average Ergodic Rate 

In this section, we derive the average ergodic rate of a typical randomly located user in the K-tier 
cellular network. Using the same approach as in (lT2l) . the average ergodic rate is given as 

K 

n = Y,'^kAk, (23) 

k=l 

where Ak is the probability that a typical user is connected to the kth tier, which is given in Lemma 
1. We denote 7?.^ as the average ergodic rate of a typical user associated with A;th tier. To compute 7?.^, 
we first consider the ergodic rate of a user at a distance x from its serving kth tier BS. The link rate 
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then is averaged over the distance x (i.e. over the fcth tier). The average ergodic rate of the fcth tier is 
thus defined as 

-JZk = E, [EsiNR, [ln(l + SINR,(a;))]] (24) 

The metric means the average data rate of a randomly chosen user within the kth tier with one active 
user per cell. It also represents the average cell throughput of the kih tier with orthogonal multiple 
access and round-robin scheduling. We first derive the most general results of the average ergodic rate, 
considering thermal noise as well as per-tier BS density \j, bias factor Bj, and path loss exponent aj. 

Theorem 2. The average ergodic rate of a typical user associated with kth tier is 

n, = ^ [ / exp I - ^ - TT Ef=i x^/s. C, it) I xdtdx, (25) 

•^k Jx>oJt>0 *- '' 

where 

C,{t) = \,pf-^{Bp+Z{e' - l,a,J,)). 
Furthermore, the average ergodic rate of overall network is 

n = J2'^7r\k [ f exp { - ^ - TT Ef=i ^'/"^ C, (t) I dtxdx, (26) 

1^^^ Jx>0 Jt>0 *■ '' 

Proof: See Appendix O ■ 

Although not closed-form, this expression is efficiently computed numerically as opposed to the 

usual Monte Carlo methods that rely on repeated random sampling to compute their results. 

We now consider the special case, where ignoring thermal noise {W = 0), and applying unbiased 

association {{Bj} = 1) and equal path loss exponents ({aj} = a) for all tiers. 

Corollary 5. When {aj} = a and {Bj} = 1, the average ergodic rate of overall network and the 
average ergodic rate of a typical user associated with kth tier are 

7^({A,}, a, 1) = 7^fc({A,}, a, 1) = T , ^, / , -dt. (27) 



(28) 



If {aj} = 4, the average rates are further simplified to 

7^({A,},4,1) = 7^.({A,},4,1)= r \ r^-r, ^^- 

Jt=o 1 + V e — 1 arctan(ve' — 1) 

Proof: If {aj} = a and {Bj} = 1, from (l25l) we obtain 

7^fc({A,},a,l) = ^ / / xexp{-nj:f^,\,pf''{l + Z{e'-l,a,l))x']dxdt, (29) 

•^k Jt>0 Jx>0 *- ^ 
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where ^ = ^^^i XjPj from (|39l) . Employing a change of variables x'^ = v and J^ e ^""dv = ^ 



give 



^-"^^'■°-^'=L l + Z(e.-l.„.l) ^"- ™ 



which is independent on the index k. From (|23] ) and Xlfc=i'^fc = 1^ we obtain TZ = IZk X]fc=i'^fc 
7^fc. This gives the resuk in (gT]). When a = 4 and {Bj} = 1, from dSl), 2:(e* - 1,4,1) 



V^ Iv / (e*-i) Thp'^^ ~ A/(e* — 1) arctan(-\/e* — 1). Combining with dTT] ) gives the desired results. ■ 
In this corollary, the double integration in Theorem 2 is simplified to a single integration, and the 
integrand is especially simple to compute. The average ergodic rate, like the outage probability, is not 
affected by BS transmit power Pj, BS density \j, and the number of tiers K. This means that adding 
BSs or raising the power increases interference and desired signal power by the same amount, and they 
offset each other. Therefore, the network sum rate increases in direct proportion to the total number 
of BSs. 

B. Minimum Average User Throughput 

We assume the orthogonal transmission, where equal time (and/or frequency) slots are allocated to 
each user one after the other in a round-robin manner. The average ergodic rate of the fcth-tier user IZk 
then means the average cell throughput in the A;th tier. For the A;th tier, the average user throughput of 
a cell is given as 

TZk = -r/, (31) 

where J\fk is the average number of user per cell of the kth tier, which is given in Lemma 2. Combining 
this equation with (|7]) and (|25l) gives 

^' = -^1 I e^p{-W-^Ef=i^^/^^QW}dtxdx. (32) 

■^fc^ Jx>oJt>0 ^ ■' 

We define the minimum average user throughput as 

Q^ min TZk, (33) 

k€{l,-,K} 



which takes the minimum value among the K values given by (1311) . The metric represents the minimum 
quality of service (QoS) that the network can provide. Since it is highly dependent on the number of 
user in a cell, i.e. cell load, the minimum average user throughput well measures the effect of biasing 
on the QoS of the HCN. 
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V. Numerical Results 

A. Accuracy of Model and Analysis 

We use a path loss at 1 meter of Lq = —38.5 dB and thermal noise W = —104 dBm (i.e. 10 
MHz bandwidth) for all numerical results. We obtain the outage probability (SINR CDF) using Monte 
Carlo methods where BSs are deployed according to the given model, and a user is fixed at the origin 
for each network realization. For each spatial realization, a SINR sample is obtained by generating 
independent Rayleigh random variables for the fading channels. 

Fig. [3] compares the outage probability of the proposed PPP BS deployment, an actual tier-1 BS 
deployment, and a hexagonal grid model for tier 1 . A total of three tiers are modeled with the lower 
two tiers (e.g. pico and femto BSs) modeled according to a PPP. We observe that the tier-1 PPP model 
is nearly as accurate as a hexagonal grid model, where the grid model provides lower bound with a 
gap less than 1 dB from actual BS deployment and the PPP model gives upper bound with less than a 
1.5 dB gap from an actual BS deployment. Similar results are also observed for two-tier case in [|25l . 
The analytic curves given from (fT6l) are remarkably close to the simulated curves for all considered 
SINR threshold, which is an advantage over [|25l that provides an exact expression for r > 1. 

B. Effect of BS Density, Path Loss Exponent, and Biasing 

We obtained numerical results of outage probability (in Theorem 1), average ergodic rate (in Theorem 
2), minimum average user throughput (from (|32l) and (|33l) ) with respect to main network parameters; 
bias factor, BS density, and path loss exponent. Although all results given in this section are for a two- 
tier HCN (e.g. macrocell and picocell), they can be applied, without loss of generality, for a general 
i^-tier HCN. Biasing effect is investigated by adjusting the bias factor of the picocell with no biasing 
of the macrocell. Since we assume the transmit power of the macro BS 20 dB larger than pico BS, 
the 20 dB picocell bias factor means the cell association where the user connects to the BS with the 
lowest path loss (or to the nearest BS for the same path loss exponents). 

Outage Probability and Average Ergodic Rate for BS density. In Fig. |4l although considering 
noise, we observe no changes in outage and rate for adding BS with different power, no biasing, 
and same path loss exponents. This means that inter-BS interference is still dominant in our two-tier 
scenario considering a typical value of BS density and BS transmit power, and Corollaries 2 and 5 
hold for the HCN. When the low-tier (pico or femto) BSs experience higher path loss, the outage 
and average rate improves as the BSs are added. This is even more optimistic than the result for the 
same path loss exponents. Intuitively, higher path loss reduces the interference between picocell and 



15 

macrocell so the picocell is more isolated from the macrocell network. The results indicate that if 
given the choice, new BSs are better deployed in an area with higher path loss. 

Outage probability and Average Ergodic Rate for Biasing. Fig. [5] shows the effect of bias factor 
and BS density. We apply the same path loss exponent for all tiers. Unbiased association results in 
the same outage/rate for all BS density as provided in Corollaries 2 and 5. As the picocell bias factor 
increases, more macro users with low SINR are associated with the picocell, which improves the 
outage and rate of picocell, but degrades those of the macrocell. In terms of the outage and rate of the 
overall network, unbiased association always outperforms biasing. Intuitively, in biased association, 
some users are associated with the BS not offering the strongest received signal, which reduces the 
SINR of the users. The results are from the condition that every cell is fully traffic-loaded (full queues 
at all times). Note that for a lightly-loaded HCN, biasing can improve the rate over the whole network 
even if it reduces the SINR of the users, since their share of the total resources (typically, time and 
frequency slots) will increase. 

In Fig. [51 for a given bias factor, deploying more pico BSs enhances the outage and rate of macrocell 
(which means a decrease in the outage probability and an increase in the average ergodic rate), because 
more macro users with low SINR and at cell edge become associated with the picocell. The outage 
and rate of picocell are also improved, since adding more pico BSs tends to reduce each picocells 
coverage area, despite the biasing. Interestingly, the BS density only slightly changes the outage of 
overall network. Although adding pico BSs increases the interference to macrocells, it also decreases 
their association probability. These two competing effects more or less cancel. Similarly, the rate is 
also not strongly affected by the BS density. 

Minimum Average User Throughput for Biasing. In Fig. [6l as the bias-factor increases, the 
minimum user rate increases at first but then decreases for a sufficiently large bias factor. Although 
the macrocell load decreases at the cost of an increase in picocell load, the average user throughput 
of picocell is still higher than that of macrocell for a smaller bias factor. For a sufficiently large bias 
factor, the average user throughput in the picocell is lower than that of the macrocell due to a massive 
number of connections to the picocell. Again, we wish to emphasize that these observed trends on 
biasing can be considered preliminary, and are heavily dependent on the cell loading. Ideally, it seems 
that a network should dynamically push users onto lightly loaded cells as a function of the current 
network conditions, and that in general, small cells will be more lightly loaded unless proactive biasing 
is introduced. 
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VI. Conclusions 

This paper developed a new analytical framework for evaluating outage probability and spectral 
efficiency in HCNs with flexibly cell association, also known as biasing. It is interesting to observe 
that the number of tiers and density of base stations at most weakly affects the outage probability 
and the average ergodic rate, and under certain assumptions, does not affect them. This implies that 
even randomly adding pico and femtocells to a network for capacity improvement need not decrease 
the quality of the network, as is commonly feared. Assuming full queues at all base stations, biasing 
deteriorates the outage and rate of the overall network by lowering the SINR, but further work on this 
topic is needed, perhaps with the assistance of the model and results developed in this paper. 
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Appendix 

A. Proof of Lemma 1 

Denote n as an index of tier associating the typical user. When Pr,fc > -Prj for all j G {1 
k, a typical user is associated with the fcth tier, i.e. n = k. Therefore, 

Ak=r[n = k] 



K},J^ 



E 



Rk 



P 



-Pr,fc(-Rfc) > maxPrj 



E 



Rk 



K 



(a) 



E 



Rk 



■ K 

n p 



Pk Bk ' ^^ 



Y\.j=h:i^k^ 



R,>{P,B,y/"^r''^ /^,(r)dr 



(34) 
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where (a) is given using ([T]). '^[Rj > {PjBjY^°'^r^^] and the PDF of Rk are derived using the null 
probability of a 2-D Poisson process with density A in an area A, which is exp(— AA). 

K K 

JJ P \Rj > (P,-5j)^/"^r"^j = Yl P[No BS closer than {PjBjfl''^r^^ in the j th tier] 

= f[ e-^>^Am?^''^r-^''\ (35) 

and 



/ij, (r) = 1 '-^^ '- = e-^^'^' 27rXkr. (36) 



^3^^^.e---^2.A 
Combining ((341), dS]), and ([361), we obtain 

Ak = 27rAfc y" r exp j-vr Ef=i,,yfc \,{P,B^f'^^r^/''^ - vrAfer^} dr. (37) 



Since Pj = 1, 5^- = 1, and aj = 1 for j = k, we obtain Eli y^fc ^i(^i^i)^^"'^^^"' + ^fc^^ 



E^i \j{PjBjY/''^r'^/''K From dST]), we thus obtain the desired resuk in ©. 
If {uj} = a, then we obtain 5?^ = 1 for all j G {1, ■ ■ ■ , K}. This gives 

Ak = 27rAfc f r exp j-vr Efli Xj{PjB^f^''r^} dr. 

Employing the change of variables P = t gives 

^k 



(38) 



Ak 






Ef=iA,(P,i?,)^/" 



Ef=i,,v.A.(^.^.)'/" + AA 



(39) 



_ „,,., ,, _ t,^ FXQ 



By applying Pj = -p- and Bj = -^ to (|39| ). we obtain the simple result in dH). 



B. Proof of Theorem 1 

From (fT3l) . the outage probability of the /cth tier is given as 

/■oo 

Cfc = 1 - / P[SINRfc(x) > r]/x,(x)da; 

9 \ /"oo 

= 1 — ^ / P[SINRfc(x) > r]xexp { -vr Ef=i \j{PjBjfl"^x'^/^^ ] dx, (40) 

where fxA'^) is given in Lemma 2. The user SINR in (fT4l) is rewritten as 7(0;) = „^p°i , where 
Q = Ei=i -^j + ^/Lq The CCDF of the user SINR at distance x from its associated BS in kth tier is 



18 



given as 



k{x) > r\ 



'k,0 



> x'^'P^^tQ] 



POO 

/ exp{-x"^p-Vg}/Q(g)dg 
Jo 

EQ[exp{-a;"'=P,-Vg}] 



.^\ 



K 






(41) 



where SNR 



W 



, and the Laplace transform of Ij is 



£,^,(x"'=P,7V) = E,Je-^'"'=^'^"^^^-] 



E*, 



exp < —X 



Qfc 



Pj'^ / , '^jA^h' 



ie^i 



(a) 



exp <^ -27rA, / 1 - £;,^.(x"'=P,r|/-"0 Z/dy 



(&) 



exp < — 27rA,- 



exp < — 27rA,- 



oo 



1- 



1 



1 + x"'=P,r|/-" 
2/ 



ydy 



dy 



(42) 



1 + {x'^kp.r)-^y'= 

where (a) is given from lfT6ll . and (b) follows because interference fading power hj ~ exp(l). The 
integration limits are from Zj to oo since the closest interferer in jth tier is at least at a distance 

Zj = (PjBjY/°'^x°'i . Employing a change of variables u = (x"'=P,r)"^/°^y^ results in 



Ci^ (x^'^Pfc-V) = exp |-7rA,pf "^Z(r, a,-, P.^x'/"^ } 



(43) 



where 



Z{T,aj,Bj) = r' 



(B. 



,l^fl^, 1 + M"j/2 



du 






2^1 



l,l--;2-^;-^ 



a, ' B, 



for aj > 2 



Here, 2-^1 [■] denotes the Gauss hypergeometric function. Plugging (l43l) into (|4TI) gives 



K 



P[7fc(x) > r] = exp 



SNR 



7r5^A,P,^/"^Z(r,a,,P,)x2/^^ 



(44) 



(45) 



Combining (|40|) and (|45l) . we obtain the per-tier outage probability in (fT5l) . Furthermore, plugging (fT5l) 
into (fT2l) gives the network outage probability in (fT6l) . 
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C. Proof of Theorem 2 

From (124)) . the average ergodic rate of the kih tier is 



nk= / EsiNR, [ln(l + SINRfc(a;))]/x,(x)dx 
'o 



27rAfc ''°° 



^ f EsiNR, [ln(l + SINRfe(a;))] x exp {-tt J^f^, \{PjBjf'^^x^'^^ } dx, (46) 
where fxA^) is given in Lemma 2. Since E[X] = Jg°° P[X > x]dx for X > 0, we obtain 

/•oo 

EsiNR, [ln(l + SINRfc(x))] = / P [ln(l + SINRfe(x)) > t] dt 

Jo 

POO 

= / P[(?fc,o>a;"^Pfe"'Q(e*-l)]dt 

e^'^Y[Ci^{x'"'Pj;\e' - l))dt, (47) 

where (a) follows from plugging r = e* — 1 in (|4T]) . From (l43l) . we obtain 

£,^.(x"^P,-i(e* - 1)) = exp |-7rA,pf "^Z(e* - 1, a„ 5,)a;2/^^ } , (48) 



(a) / ^sLzl 



with 



Z(e* - 1, a,-, R) = (e* - 1) "^ / rrrfw 



2(e* - 1)r'/"^-' 
a,- — 2 



2 . o 2 . l-e' 



l,l-^;2-^ 



Qj "J £>j 



Plugging (gSj) into (147]) gives 



ic 



(49) 



EsiNR, [ln(l + SINRfc(x))] = J exp J -^^ - tt ^ A,i52/"^x2/"^Z(e* - 1, aj, Bj) [ dt. (50) 

Combining (|46l ) and (|50l ). we obtain the average ergodic rate of the fcth tier in (|25] ). Furthermore, 
plugging (l25l) into (l23l) gives the average ergodic rate of entire network in (|26l) . 
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Received power from pico BS 
+ Bias factor (>0) 



Received power 
from macro BS 




Received power from pico BS 
> Received power from macro BS 



Received power from pico BS + Bias factor 
> Received power from macro BS 



Fig. 1. Concept of flexible cell association with positive bias factor in a two-tier macro/pico HCN. 
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Fig. 2. Example of downlink HCNs with three tiers of BSs: a high-power macrocell BS (red square) is overlaid with successively 
denser and lower power picocells (blue triangle) and femtocells (magenta circle). Both figures have the same BS position. Black lines 
show the border of cell coverage. The pico/femtocell with positive biasing (lower figure) gives a larger coverage than no biasing (upper 
figure). 
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Fig. 3. Outage probability comparison in a tliree-tier HCN ( A2 
{3.8,3.5,4}, and {Pi, Pa, P3} = {53, 33, 23} dBm, ). 



2Ai,A3 = 2OA1, {Bi,B2,B3} = {1,1,1}, {01,02,03} 
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Fig. 4. Outage probability and average ergodic rate for varying density of picocells in a two-tier HCN {K = 2, Ai 
{Bi,B2} = {1,1}, and {Pi,P2} = {53,33} dBm). 
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Fig. 5. Outage probability and average ergodic rate for varying bias factor of picocells in a two-tier HCN (K — 2, Bi = 1, Ai — J^^ , 
{01,02} = {3.5,3.5}, and {Pi,P2} = {53,33} dBm). 
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Fig. 6. Minimum average user throughput for varying bias-factor and density of picocells in a two-tier HCN {K — 2, Bi = 1, 
^1 = dw^' {^1.^2} = {53,33} dBm). 



